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Nonlinear elliptic-parabolic problems 

Inwon C Kini*and Norbert Pozar''' 



Abstract 



o 

» I We introduce a notion of viscosity solutions for a general class of elliptic-parabolic 

d phase transition problems. These include the Richards equation, which is a classical 

^^ model in filtration theory. Existence and uniqueness results are proved via the com- 

parison principle. In particular, we show existence and stability properties of maximal 

^^ and minimal viscosity solutions for a general class of initial data. These results are 

^~~^ new even in the linear case, where we also show that viscosity solutions coincide with 



the regular weak solutions introduced in AL 

"*^. 1 Introduction 



Let Q C M" be a smooth bounded domain and T > 0. Let us denote Q ~ ft x (0,r]. We 
are interested in the following problem: find a function u{x,t), u : Q ^>- R, that solves 

'dtb{u)~ F{D'^u,Du,u)^0 in g, 
^ u = g=-l ondnx[0,T], (LI) 

■^ [m(-,0)=mo on ft, 

(N 

CN where Du denotes the spatial gradient of u, D^u is the spatial Hessian, and F{M,p,z) : 

^^ Sn X K" X M — > M is a fully nonlinear, uniformly elliptic operator (see Section [lT] for precise 

C^ assumptions on F) . For the function 6 : M ^- ]R we assume that 

O 

OQ (a) b is increasing and Lipschitz, 

L; (b) b{s) = for fo < 0, 6 e C(M) n C\[0,oo)), 

k> (c) there exists a constant c > such that b'{s) > c for s E (0, oo). 

5—1 

^ The nonlinear operators F we consider include: 

• the uniformly elliptic operator of non-divergence form 

F{D^u, Du, u) = - tT{A{Du)D^u) + H{u, Du), (1.2) 

where A satisfies the uniform ellipticity condition 

A|gP < A{p)q ■ q < A\q\^ for all p,qe M", for some A, A > 0, (1.3) 
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and as well as the Bellman-Issacs operators arising from stochatic optimal control and 
differential games 

F{D^u,Du,u) ^ inf SMp{C"^u}, 



where C"^ is a two-parameter family of operators of the form (1.2) satisfying (1.31 
we refer to [CIL FS| for further examples. 



• a divergence-form operator; to simplify our discussion, we restrict our attention to 
operators of the form 

F{D^u,Du,u)=V ■{■9{b{u))Du), (1.4) 

where ^ e C^([0, oo)) is a positive function. The class of operators given in ( |1.4[ ) is 
of particular interest since in that case the problem (1.1) represents the well-known 



Richards equation, which serves as a basic model for the filtration of water in unsatu- 
rated soils (see e.g. [DS^MR^Rj). 

Our aim is to study the well-posedness of (1.1[). Note that, due to the regularity theory 



for uniformly elliptic nonlinear operators ( CC 



W2 



), solutions of (1.1) satisfy interior C^'" 



estimates in the sets {u > 0} and {u < 0}. Hence the challenge in the study of the problem 
lies in the behavior of a solution near the transition boundary between the positive and the 
negative phases: for example, as illustrated in Example |1.2| discontinuities of solutions in 
time across the set {u = 0} are generic. 

Problem ( |1.1[ ) can be understood as the limiting equation for the evolution of two phases 
with different time scales of diffusion and with the threshold value at u = 0. In particular, as 
verified in section [6j the problem can be viewed as a singular limit of a family of uniformly 
parabolic problems ( |6.1[ ), where b{s) in (1.1) is regularized. Hence it is expected that a 
maximum principle holds for the solutions, and the theory of viscosity solutions may be 
applicable for the study of pointwise behavior of solutions near the transition boundary 
{u = 0}. This is indeed our approach: in this paper we will introduce the notion of viscosity 
solutions for ( |1.1[ ) and discuss existence, uniqueness and stability properties, and compare 
them to the notion of weak solutions (see the discussion below). Such results have been 



established for Stefan- type problems (see CS KP , for example), but significant challenges 



in the analysis arise due to the implicit nature of the boundary motion law in ( 1.1 1 and the 



generic nature of discontinuities in u (as in Example 1.2 ); see an extended discussion on this 
in section [3] We aim to present the proof of the comparison principle for fully nonlinear 
operators in more detail, both to illustrate the flexibility of the viscosity solution approach 
and to make the results readily available for applications in a general context. 

It should be pointed out that our approach treats the transition boundary T := d{u > 0} 
between the elliptic and parabolic regions as a "free boundary" and constructs barriers based 
on its movement. Hence our approach may not be optimal for other problems such as fast 
diffusion (when 6'(0) = 0) where the transition boundary moves with infinite speed (see [V] , 
for example), but strong regularity properties for the solutions are expected (see [BH|). 



Although our presentation is mainly focused on the problem with fully nonlinear oper- 
ators whose structure assumption requires a linear growth in terms of \Du\, the results of 
this paper extend to the case of divergence- form operators of the form (1.4 1 as well. We 
point out the necessary modifications in the text where appropriate. 



Literature review 



Let us briefly discuss previously known results on (1.1 1: all of them concern divergence- 



form operators including (1.4 1. The weak solutions are defined via integration by parts in 
the important paper of Alt & Luckhaus |AL , which shows existence ( [AL[ Theorem 1.7]) 
of weak solutions for the general class of elliptic-parabolic phase transition problems with 
divergence-form operators. Uniqueness results are, however, rather limited. For F given 



by (1.4), when ^(s) is a positive constant (that is, when F is linear), the authors prove 
that regular weak solutions (in the sense dtb{u) g L^) can be constructed by the Galerkin 
method ( AL Theorem 2.3]), and show that regular weak solutions satisfy the comparison 
principle and are unique given the initial data (lALl Theorem 2.2]). Indeed, when ^(s) is 



a positive constant, it is known that b{u) is continuous in a local setting as long as u is 
bounded (see DG| ). Continuity of b{u) seems to be the optimal result for this problem with 
general initial data, since u may become discontinuous in time in the elliptic phase (see 



Example 1.2, and also [DG| ). The proof of continuity in DG is based on the weak Harnack 
inequality, making use of the linearity of the elliptic operator with respect to u. We also 
refer to 



BW 



and ICl, who use an entropy solution approach to define weak solutions as 
well as to prove comparison principle in L^ for the relevant model; this approach, while 
powerful, does not fit into our setting where we have a non-vanishing elliptic phase. 



Even for the quasi-linear F given in (1.4), there are no uniqueness or stability results 



except for the aforementioned linear case and for the one-dimensional case (see BH vDP ): 
this serves as a motivation of our analysis in this paper. In one dimension, Mannucci 
& Vazquez MV studied viscosity solutions of (1.1) for divergence-form operators. Their 
approach avoids possible complications at the transition boundary d{u > 0} by relying on 
previously known regularity properties of weak solutions in one dimension. 



Summary of the main results 

In this section we summarize the main results obtained in this paper. In all statements F 
is assumed to be either a fully nonlinear operator satisfying the assumptions in Section |1.1| 
or a quasilinear divergence- form operator of the form (1.4). 



Our main theorem is the following comparison principle: 



Theorem 1 (Theorem 3.1 and Theorem 3.24). Let u be a viscosity subsolution and v a 
viscosity supersolution of ( 1.1 ) on Q = Q x (O^T] for some T > 0. If u < v on the parabolic 
boundary dpQ, then u < v on Q. 

Equipped with the comparison principle, we use Perron's method to show the following 
existence and stability theorem. 



Theorem 2. For initial data uq E V (see the definition of V in section^ by (4.1 )-(4.2)J, 
the following holds: 



(a) (Theorem 4-S) There exists a minimal and a maximal viscosity solution u and u of 
(1.1) with initial data Uq. 



(b) (Theorem 4-3) u and u are stable under perturbations of initial data with appropriate 
ordering. 



(c) (Corollary 6.3) u and u can be obtained as a limit of solutions solving the regularized 



parabolic equation (6.1) 



Aforementioned theorem states that the maximal and minimal viscosity solutions are 
stable. Unfortunately, we are only able to show the uniqueness of general viscosity solu- 
tions (i.e. the coincidence of minimal and maximal viscosity solutions) in several restricted 
settings. The coincidence of the minimal and the maximal viscosity solutions with general 
initial data remains open, except for the linear case. 

Theorem 3. For given initial data uq Cz V, the following holds: 



(a) (Theorem 5.6) If F is linear, i.e. if F{M,p, z) = F{M) = tv M , then there exists a 



unique viscosity solution u with initial data uq, and u coincides with the unique weak 



solution defined in AL 



(h) (Theorem 5.8) If uq is either star-shaped or if u decreases att = 0, then there exists 



a unique viscosity solution of (1.1). 



1.1 Remark. As mentioned above, our approach may not be optimal if b{s) degenerates at s = 0. 
On the other hand, we expect that our approach can be extended to non-Lipschitz b{s) and produce 
results similar to the above. The difficulty in the analysis lies in the corresponding degeneracy of 
the elliptic operator in the positive phase, when we write the equation in terms of b{s). 

1.2 Example (Discontinuous solution). Here we briefly discuss an example which illustrates dis- 
continuities in the solutions. Set 

b{u) = u+ :— max(0,u), 

and consider ( |1.1[ | with negative boundary data, and initial data that are positive on some open set. 
As the solution evolves, the positive phase disappears in finite time, and then the solution jumps 
to the stationary solution. Nevertheless, one expects b{u) to be continuous. We refer the reader to 
lALl p. 312] for an explicit formula. 



1.1 Assumptions on the nonlinear operator F 

Let Sn be the space of symmetric n x n matrices. For given < A < A, we define the Pucci 
extremal operators A4^ : 5„ -^ M as in CC W 



M + {M) = sup tr AM, M^{M) = inf trAM, (1.5) 

Ae[A/,A/] Ae[xiM] 

where [A/, A/] = {A ^ Sn : \I < A < A/}. Alternatively, the Pucci operators can be ex- 
pressed using the eigenvalues ei, . . . , e„ of matrix M: 

M+{M) = A ^ e, + A ^ e„ M' {M) = A ^ e^ -f A ^ e,. 

ei>0 ei<0 ei>0 ei<0 

With the Pucci operators at hand, we shall assume the following structural condition on 
the operator 

F{M,p, z) : 5„ X M" X M ^ M : 

(i) There exist < A < A and Sq.Si > such that 

M-{M -N)-S,\p-q\- So \z-w\< F{M,p, z) - F{N, g, w) 

<M+{M-N)+6i\p-q\+5o\z-w\ 

for all M, N eSn,p,qe M" and z,w eR. 



(ii) F is proper, i.e. 

z i-^ F{M,p,z) is nonincreasing in z. (I-''') 

(iii) Finally, to guarantee that u = is a solution of both the parabolic and the elliptic 
problems, we assume that 

F(0,0,0) = 0. (1.8) 

1.2 Notation 

• In this paper, we work in a fixed space dimension n > 2. For a point x G M" and time 
i e M, the pair {x, t) E M" x M represents a point in space-time. 

• For given r > 0, we define the open balls 

Br := [{x,t) : \x\^ + \t\^ < r^}, B," := {x : \x\ < r}, 

the space disk 

Dr := Bl! X {0} = {{x,Q) : \x\ < r}, 
and the flattened set 

Er:^{ix,t):\x\' + \tf<r'}. 

Finally, we define the domain that is used in the definition of regularizations of solutions. 

Here + is the Minkowski sum. Note that S^ € C^ (in contrast with Dr + Br^ which is only 
Ci'i). 

• We ask the reader to forgive a slight abuse of notation: 

dDr ■.= Wr-Dr = dB^! X {0} = {(a;,0) : \x\ = r}. 

• It will also be advantageous to introduce the (open) top and bottom flat pieces of d'E.r, 

d'^Er := {{x, r) : \x\ < r}, d±Er = {{x, -r) : \x\ < r}, 

and the (open) lateral boundary of 5^, 



dL^r '■= d^r \ 9^S,. U d±^r- 

These sets are sketched in Figure [l] 

• The translation of a set yl e M" x M by a vector (x, t) e M" x M will be denoted as 

A{x,t) := {x,t) + A. 

The translation A{x) of a set A C M" is defined similarly. 



tt 




d'^-r i\x\- 



Dr 



d±^r 




t = r 



t = -r 



Figure 1: The boundary dEr of the set 5,. 

• We will often need to consider timeslices (cross sections at a fixed time) of a given set. To 
simplify the notation, let us define the timeslice of a set A C M" x M at time t as 

A\^ := {x : {x,t) G A}. 

We often write Af if there is no ambiguity. 

• Let i? C M" X M. Then USC{E) and LSC{E) are respectively the sets of all upper semi- 
continuous and lower semi-continuous functions on E. For a locally bounded function u on 
E we define the semi-continuous envelopes 



inf V, 

veusc(E) 



U*.E 



sup V. 

veLSC(E) 



(1.9) 



These envelopes are used throughout most of the article with E = Q, and therefore we 
simply write u* and u^, if the set E is understood from the context. 



2 Viscosity solutions 

In this section we define the notion of viscosity solutions of problem ( |1.1[ ). Formally, viscosity 
solutions are the functions that satisfy a local comparison principle on parabolic neighbor- 
hoods with barriers which are the classical solutions of the problem. We refer the reader to 
CIL CV CS KP| and references therein for other examples of this approach. 



2.1 Definition (Parabolic neighborhood and boundary). A nonempty set E C K" xM. is 
called a parabolic neighborhood if E = U H {t < t} for some open set U C M" x K and 



some T g M. Let us denote dpE := E\E, the parabolic boundary of E. (See 
more general definition.) 



W 



for a 



2.2 Definition (Classical subsolution) . Let E be a parabolic neighborhood. Function (p is 

trhood E if ip G C{U) on 
and the following holds: 



called a classical subsolution of problem, ( 1.1 1 in a parabolic neighborhood E if ip G C{U) on 
an open set [/ G M" 



< M such that E = Ur\{t < t} for some t € 
(z) ^ e CljiJ^^O}) and CljiW^^), 
(a) {ip = 0} c d{ip > 0} n d{<f < 0} and \Dip^\ > on {ip 



0}, 
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(m) b{if)t - F{D'^ip, Dip, ip) <Q on {if > 0} and {ip < 0}, 
(iv) \Dip+\ > \Dip-\ on {ip = 0}. 
Here {ip > 0} := {(.T,t) G U : ip{x,t) > 0} etc., and 

D^^{i,T):^ lim D^{x,t). 

(x,t)-^(C,r) 
(a:,t)e{±¥'>0} 

We say that ip is a strict classical subsolution if the inequalities in (Hi) and (iv) are 
strict. 



Classical supersolutions are defined similarly by flipping the inequalities in Definition |2.2| 
(iii)-(iv). 

At last, we define viscosity solutions. Note that we set g = — 1 on dVl x [0, T] throughout 
the paper. 

2.3 Definition (Viscosity subsolution). Let Q = Vtx{Q,T] he a parabolic cylinder. Function 



u £ USC{Q) is a viscosity subsolution of (1.1 1 in Q if u{-,0) < uq on Q, u < g on 



dfl X [0,r], and if u < ip on E for any strict classical supersolution ip on any parabolic 
neighborhood E G Q for which u < ip on dpE. 

One can define viscosity supersolutions accordingly, as a function in LSC{Q), by switch- 
ing the direction of inequality signs in the previous definition. 

2.4 Definition (Viscosity solution). Locally bounded function u is a viscosity solution of 



(1.1) on Q ifu*'^ is a viscosity subsolution on Q and u^q is a viscosity supersolution on 



Q. 

2.5 Remark. We only test viscosity solutions by strict classical barriers. It is therefore possible 
to narrow the choice of E in Definition |2.3| to only include parabolic cylinders of the form Q' = 
n' X (ti, ^2] C Q, where fi' has a smooth boundary, instead of all parabolic neightborhoods. Indeed, 
suppose that _E C Q is a parabolic neighborhood, (^ is a strict classical supersolution on E, u < ip 
on dpE, but u > ip at some point in E. Define r := sup {a : u < p on E f] {t < a}} G R. The set 

A:= {x : {x, t) £ E,u> p} 

is compact and therefore 5 := dist(A x {t},9p-B) > 0. Define the parabolic cylinder Q' = [A -\- 
Bs/2) X (r — (5/2, r]. Clearly Q' G E and u < p on dpQ' . The boundary of A + 5^/2 can be easily 
regularized. This observation will be useful in section [5] where we will show that regular weak 
solutions are viscosity solutions. 

3 Comparison principle 

This section is devoted to the proof of the following "weak" comparison principle. 



3.1 Theorem. Let u be a viscosity subsolution and v a viscosity supersolution of (1.1 1 on 
Q ^ CI X (0, T] for some T > 0, and assume that u < v on dpQ. Then u < v on Q. 

To simplify the exposition of the proof of this theorem, we shall assume that T ^ 00. In 
fact, it is always possible to extend u and v from Ct x (0,T] to fi x (0,oo). Moreover, we 



will only consider b(u) of the form b{u) = u+ := max(w, 0). The problem (1.1) with a fully 
nonlinear operator F and a more general b can be always rewritten in this way. Indeed, the 
differentiation b(u)t — b' {u)ut for u > is justified by the regularity of b, and h' iu) can be 
absorbed into F. 



Heuristic arguments 



The comparison principle for classical subsolutions and supersolutions u and v of (1.1) can 
be proved using the following formal argument: we would like to show that 

{u{-,t) > 0} C {u(-, t) > 0} for aU t > 0, (3.1) 

since then the conclusion follows due to the standard elliptic and parabolic comparison 



principle. Hence suppose (3.1) fails at some time. Since \Du\, \Dv\ > on the boundary 
of their respective positive phases, the sets {u{-,t) > 0} and {v{-,t) > 0} have smooth 
boundaries and evolve continuously in time with respect to the HausdorfF distance. Therefore 
it follows that there exists the first time io > where d{u{-, t) > 0} intersects d{v{-,t) > 0}, 
let's say at x = xq. Since 

{w > 0} n {i < to} C {i; > 0} n {i < to} and uq < ^o, 

the comparison principle applied to u and v respectively in the sets {u > 0} and in {v < 0} 
yields that u < v up to t — tg. In particular, since u(xo, to) = i^i^a, to) — 0, this means that 

\Du+\{xo,to) < \Dv+\ixo,to) and \Dv-\{xo,to) < \Du~{xo,to)\. 

Furthermore, due to the regularity of d{u{-, t > 0} and d{v{-, t) > 0} and Hopf's lemma for 
uniformly parabolic and elliptic operators, it turns out that the above inequalities are in fact 
strict. This contradicts the flux-matching condition for the classical sub- and supersolution 
(Definition |2.2[ iv)), i.e. the fact that 

\Du+\ > \Du-\ and \Dv+\ < \Dv-\ at (xo,to). 

Unfortunately, the rigorous version of the above heuristic argument is rather lengthy. 
There are many difficulties one faces in the general setting, where u and v are merely 
semi-continuous functions. As it is always the case in the theory of weak solutions, one 
should translate the above heuristics onto appropriate test functions or, to be more precise 
in our case, barriers. Our argument relies on a certain regularization procedure (see sub- 
section 3.2) which ensures that, at a contact point (xo,to) of the regularized solutions, the 
phase boundaries of each regularized solution are both locally C^'^ in space. Such regularity 
of the phase boundary would enable us to construct appropriate barriers which would allow 
one mimic the heuristic argument above. This technique was pioneered by Caffarelli and 
Vazquez (|CV|) in their treatment of viscosity solutions for the porous medium equation. 



It was later applied to several one phase free boundary problems in BV K] K2| , and later 



extended to two-phase Stefan problems by the authors in [KP| (see also [CS 



In contrast to the aforementioned results for free boundary problems, the analysis of our 
problem presents several new challenges. The most obvious challenge arises from the flux 



matching condition of ( 1.1 1 on the transition boundary. While the regularization procedure 
provides regularity information in space variables, one should still show the flnite propaga- 
tion property of the phase boundary. In the aforementioned free boundary problems, the 
free boundary motion law prescribes the normal velocity of the free boundary in terms of the 
gradient of the solution, which links the space regularity to the time regularity of the solu- 
tion. Here one does not have such a direct relationship between time and space regularity of 



the transition boundary. Indeed, the flux matching condition of ( 1.1 ) turns out to be more 



delicate than the prescribed gradient condition in flame propagation type problems, since 
one has to account for the possibility that the fluxes from both sides will either degenerate to 
zero or diverge to infinity. This is overcome by the observation that with the regularization 



we can rule out the scenarios of a sudden extinction of the elliptic phase (Lemma 3.12 ), or a 



sudden shrinkage/discontinuous expansion of the parabolic phase (Lemmas 3.13 and 3.19) 
at a contact point. Note, however, that even with regularized solutions, the elliptic phase 
might instantly become extinct away from a contact point. 

We point out that, to allow for the regularization procedure, the strict ordering of u and 
V on dpQ in the statement of Theorem |3.1| is necessary. We also point out that a proof via 



doubling of variables, a classical tool in the theory of viscosity solutions (see CIL ), is not 
available for phase transition-type problems, including ( |1.1[ ). 

We shall present the proof by spliting it in a number of smaller intermediate results, 
which will be later collected in 33.61 below. 



3.1 Properties of solutions 

In this section we clarify what we mean by the parabolic and the elliptic problems. We 



refer the reader to CIL for the precise definitions and the detailed overview of the standard 
viscosity theory. 

3.2 Definition. We say that w{x, t) is a solution (resp. suhsolution, super solution) of the 
parabolic problem in an open set Q' C M" x R if w is the standard viscosity solution (resp. 
suhsolution, super solution) of 

Wt-F{D^w,Dw,w)^0 mQ'. 

Similarly, we call w{x) a solution (resp. sub /super solution) of the elliptic problem in an 
open set Q,' C K" if w is the standard viscosity solution (resp. sub/supersolution) of 

-F{D'^w, Dw, w)^0 in ft'. 

Now with the help of the previous definition, we can rigorously express the intuitive fact 



that the problem (1.1) is parabolic in the positive phase and elliptic in the negative phase. 



3.3 Lemma. If u is a viscosity subsolution of (1.1) in Q then u+ is a suhsolution of 
the parabolic problem in Q. Similarly, if v is a viscosity supersolution of ( |l.l[ ) in Q then 
— u_(-, t) is a supersolution of the elliptic problem in fi for each i > 0. 

Proof. 1. The first claim follows easily since u+ = max{u, 0} and solves 

ut- F{D^u,Du,u) = 



(see (1.81). 



2. To prove the second claim for — u_ := min{w,0}, suppose it is not a supersolution 
of —F{D^u,Du,u) = at time t = to > 0. Let us denote C(a;) = —v^{x,to). Then 
there is a function (p E C'^, f < 0, and an open (space) ball B C {C < 0} such that 
—F{D'^tp, Dip, ip) < on B, (f < C, in B and ip <( on dB, but (y3 = C at some point in B. 



Due to the continuity of F, there exists rjo > such that ~F{D'^if, Dip, (p — rj) < m B for 
all 77 e [0, 770]. For any S > we set 

Qs := B X {to-S,to]. 

Since v is lower semi-continuous, there exists small S > such that Qs CI {v < 0}, ip < v{-, t) 
on dB for t G [io ~ <5, to], and z;(a;,i) — ^{x) > — r/o on Q^. Let us define the barrier 
ili{x,t) = (f{x) + ft. Observe that -F{D'^'tlj,D'ip,ip) < and ip <0 in Qs and therefore it 



is a strict classical subsolution of (1.1). Furthermore, iJj < v on dpQs, while ip = v at some 
point in Qg. This contradicts the fact that u is a supersolution of (1.1 1. D 



3.2 Regularizations 

In this section, we define regularizations of the subsolution u and the supersolution u, and 
prove some of their properties that are applied in the proof of the comparison theorem. 
For given r £ (0, 1), define the regularizations 



(3.2) 



Z{x, t) :— sup u, 
W{x,t) :=_inf V. 

SAx,t) 

Functions Z and W are well-defined on the parabolic cylinder Q^, 

Qr := {{x,t) e Q I 'E.r{x,t) C Q} = r^r X (r, 00), 
where 

nr-={xen\ dist(a;, Q") > r + r^/^j. 

3.4 Remark. Jlr has the uniform exterior ball property (with radius r + r^'^) and therefore it is 
a regular domain for the elliptic problem. 

3.5 Remark. The main advantage of regularizing over the set H,- instead of the ball Br is that the 
appropriate level sets of Z and W have both space-time and space interior balls, see Proposition |3.10| 

3.6 Remark. Defining H = D + E, instead of D + _B as in the previous paper 



consequence that the parabolic boundary of H is not C^'j at the top flat piece 9^S, and thus 



1,1/2 



KP , has the 



it is not a regular set for the heat equation (see CS ). We can therefore expect that the gradient 



of a positive caloric function in ^ which is zero on the l ateral boundary of E will blow up at the 



lateral boundary as we approach 3^S, see Lemma 3.13 A similar effect of vanishing gradient is 
expected when the solution is positive on H'^ but zero in H, see Lemma [3.16| This is the necessary 
new ingredient required in the proof of the "finite speed of propagation" in Lemma [3.19[ 

3.7 Proposition. Suppose that v is a vise, supersolution on Q ~ fl x (0,oo). Then W is 
a vise, supersolution on Qr- Similarly, if u is a vise, subsolution on Q then Z is a vise, 
subsolution on Qr- 

The strict separation of u and v on the parabolic boundary of Q allows us to separate 
Z and W on the parabolic boundary of Qr. 

3.8 Proposition. Suppose that u G USC and v G LSC in Q sueh that u < v on dpQ. 
Then there exists rg > such that Z < W on dpQr for all < r < rg. 
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Proof. Standard from semicontinuity. D 

Arguably the most important feature of regularizations Z and W is the interior ball 
property of their level sets, as well as of the time-slices of their level sets. We formalize this 
fact by introducing the notion of dual points. 

3.9 Definition. Let r > 0, u E USC{ft) and let Z be its sup- convolution. Let P E Q^.. We 
say that P' E Q is a sup-dual point of P with respect to u 

if P' e E{P) and u{P') = Z{P). 

Let us define n"(P) to he the set of all sup-dual points of P with respect to u, 

n"(P) = {p' e s,(P) : u{p') = z{P)). 

Similarly we can define inf-dual points Tiv{P) for v G LSC{fl) by Il^{P) :— n~"(P). 

In what follows, we shall use the following convenient notation for various level sets of u 
and Z: 

{w > 0} = {{x,t) e Q : u{x,t) > 0}, 
{m < 0} = {(x,t) e Q : u{x,t) < 0}, 

which is contrasted with 

{Z > 0} = {{x,t) e Wr ■■ Z{x,t) > 0}, 
{Z <0} = {{x,t) e Qr ■■ Z{x,t) < 0}. 

Sets {w < 0}, {w > 0}, {W < 0} and {W > 0} are defined in a similar fashion. This choice 
guarantees that sets with > and < are closed, while sets with < and > are open. 
We first make a few simple observations about 11" and 11^. 

3.10 Proposition. Let u e USC(Q). Then for all P eQ^: 

(z) n"(p)^0, 

(a) S^(P') nQ^diZ > Z{P)} for all P' £ n"(p), 

{lii) ifPe d{Z > 0}, then n"(P) C aS^(P) n d{u > 0} and S^(P) C {u < 0}. 

Since the closed sets {W < 0} and {Z > 0} have closed space and space-time interior 
balls at each point, they don't have any points that are isolated from their interior: 

3.11 Lemma. The level sets of the functions W and Z defined above have the following 
properties: 



int{M^<0} = {W^<0}, int{Z>0} = {Z>0}. 

Moreover it is true for every time-slice t (with space closure and interior). 

Proof. This follows from the interior ball property in Proposition |3.10| Fix time t and 
for simplicity define E = {W < 0}^. Pick any point x £ E. Then there is an open ball 
B — Br{y) such that x £ B and B C E. But B is open so P C intP and therefore 
X £B C intE. n 
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d{w < o}„ 




d{v<o}. 



d^r ; <9ri 



Figure 2: Space cross section of the situation in the proof of Lemma 3.12 



Another important property of the regularized supersolution W is that each point of the 
time-shce {W < 0}^ is connected to the boundary dflr with nonpositive values of M^ by a 
wide "trunk" of finite length: 

3.12 Lemma. Let W be the inf-convolution of a supersolution v. Then for every (^,(t) S 
{W < 0} there is a piecewise linear continuous curve 7 : [0, 1] — >■ fir, with finite length such 
that 7(1) G fir and 

^e U Br/2{i{s))nnrc{w<0}\,^,. 

se[0,l] 



Proof Pick P = (^,0-) e {W < 0}. By Proposition |3.10[ there exists a point P' = {C,cr') 
on n„(P) such that {^,cr) £ Er{^',a') C {W < 0}. Therefore we can also find ^ such that 
£. G dBr/2{i) and 5^/2(0 C {W^ < 0}^. The situation is depicted in Figured 

Let us denote H ^ {v < 0}|j^^,. Since H is compact we can select a finite subcover 
from the open cover 

i/ C U Br,2{x). 

xeH 

Let us denote the balls in the finite subcover by B^ ^ . . . , B^ , with centers xi, . . . , x^. 
Suppose that there exists a permutation (ji, . . . ,jk) of (1, ... , k) and (7 G N such that 

dfl\J [j B^' I = 0. 

q<l<k J 

Let us denote G = Ui</<(7 ^"" and C = HDG. Note that G is open and C is nonempty. By 
definition, u > on dGx {a'} and therefore there is (5 > such that w > on dGx {a' — 6, a'] 
{{v > 0} is open by u G LSC). In particular, since — w_(-,t) is a supersolution of the 




elliptic problem for every t (Lemma 3.3 1, an application of the elliptic comparison principle 
(Proposition |A.l ) yields that v > in G x (u' — (5, a'] . And a straightforward barrier 
argument (barrier constructed in Lemma B.4) shows that in fact u > on G x {cr'}. This 



is a contradiction with v{xj^,a') < at Xj^ G G. 
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Figure 3: Iterations in the proof of lower bound in Lemma 3.13 



Therefore we conclude that we can find ii, . . . ,im distict with ^' e i3'\ -B*'" fl d^ 7^ 



and 



B*' n S*'+i 7^ 0, 1 < ? < m. 



Since we observe that Br{xi) C {VF < 0}^ for all 1 < i < A:, we can choose the curve 7 as 
the piecewise linear curve connecting points ^, a;^^ , . . . , Xi^ . D 



Now we present two results, Lemmas |3.13| and |3.16[ that justify the shape of the domain 
Sr chosen in the definition of the regularizations Z and W . 



3.13 Lemma. Let v he a vise, supersolution of (1.1) on ^, ^ — :^r{^,a) for some (^, cr) G 



M" X R, and V > in E. Then there exists f G C([0,r]), /(O) = 0, / > on (0, r] and 
'^^^ ^- 00 as s —i' 0+ SMc/i t/iat 

f(2:,i) > fir- \x- CD for \x-C\ <r, t e [a + r/2,a + r]. 

Proof. We only prove that v{x,a- + r) > f{r — |.x — ^|). The full result follows from the 
simple observation 



^r/2{x,t) C Er for \x\< -, te 



r r 
2' 2 



Let us fix a point (C, f + r) ^ Dr{C,(J + r). Since the argument is invariant under 
translation and space rotation, we can assume that (CjCt -\- r) = (0,0). Let s = r — |C — ^|. 
Our goal is to show that ^(0, 0) = v{C,, u + r) > /(s) for some function /. The situation is 
depicted in Figure [3] 

It is straightforward to estimate the distance of the lateral boundary OlE at each time, 
9S|j, from the origin x = for t G (— r, 0): 



dist(0, 9Sy = s+ ^r2 - {t + r^ = s + ^/ -2rt - i^ > s + ^/-rt. (3.3) 

Since v G LSC is positive in S, it has a positive minimum on the compact set 

K := {{x,t) ■.\x~£,\<r+ ^/^t, t G [-7r/8, -r V256] } 
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as iiT c ^. 
Let us set 



7 := mm 



1 



and let us define the barrier 



Lp{x, t) 



16nA + 8(5i + 2(5o ' 



~t-A\x\' + l 
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(3.4) 



+ 



Clearly (^ is a viscosity subsolution of the parabolic problem due to (1.6 1, and (/3(0,7) 
^(/3(0, 0) = ^. Let us define a sequence of parabolic cylinders 



Q, = Ba^ X (/ij - (1 + 7)a', /ij), i > 0, 



where 



i=0 



0. 



flj's are defined in accordance with (3.3) as 



= ^/-rhj + s, j > 1, 



Oq = mm 



{'■{.} 



Writing ft.j+i = hj — a^, we can derive the recurrence relation 



flj+i 



^ra^j + (oj - s)3 + s, j > 0. 



(3.5) 



Note that a^-'s were chosen in such a way that Qj C ^ for all j for which /ij+i < r due to 



(3.3) 



First assume that s < r/16. From (3.5) we estimate 



a,+i > r^^^af\ j > 0, 



(ly 



which yields 

Oj > r I - j — > r as J — >■ cx). 

We want to estimate the minimal k = k{s) that guarantees 
dsQk = Da, (0, hk-{l+ j)ak) C K. 



(3.6) 



We first observe that if a^ > | + s for some j > then k < j. Therefore (3.6) yields the 
upper bound 



k< 



log 



I^. 



log I 



(3.7) 
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To show that such k indeed exists, suppose that hj > — r^/8 for some j. Then qj < r/2+s 
and 



/i,-(l + 7)a,'> -'^-2(^+5) 



' 2rs - 2s^ > - — 



since s < r/16 and 7 < 1. 

If s > r/16 then we set k = 0. 

Next we iteratively define, for j = fc, . . . , 0, the rescaled barriers 

(Pj{x,t) — Kjip{aJ^x,aJ^{t — hj) + 1 + 7), 
where 



{miiiK V > j = k, 

™nSa,/2 '/'j+i(-, hj+i - ja'j) j < k, 



and the functions tpj, j — fc, . . . ,0, are the unique solutions of the parabohc problem with 
boundary data ijjj — ipj on dpQj. We observe that ■0^ > ipj in Qj due to the parabolic 



comparison (Prop. A. 2 1 since ipj is a subsolution of the parabolic problem. 



Now the parabolic Harnack inequality (Prop. A. 4 applied with ii = 7 and ^2 = 1 + 7) 
yields 

i^j-i = „ inf i'ji-^hj - 7a2_i) > inf Vj 

> - sup ipji-^hj+i) > — sup Vj(-,/ii+i-7ai) (3.8) 



i/2 "^-Daj/2 

1 

2^' 



-Kj. 



Finally, we realize that, by definition, v > ipk — ipk on dpQk and thus an application of 
the parabolic comparison shows that v > V'fc in Qk- Since we chose ipj so that ^j < ipj+i 
on Qj n Qj+i, we can apply the comparison principle iteratively to conclude that v > ipj in 
Qj for j = fc, . . . , 1. The case j = 0, however, has to be considered separately, because we 
only know that u > in S, i.e. in Qq n {i < 0}. Therefore the parabolic comparison can 
only show that v > ipo in Qq Ci {t < 0}. Nevertheless, since Vo(0,0) > and v is bounded 
from below on S, a straightforward barrier argument using a strict classical subsolution of 



( 1.1 ), that can be constructed in the form similar to (3.4), extended in negative phase using 
a strict subsolution of the elliptic problem, shows that w(0, 0) > ■00 (0, 0). 

Now let a = ^. Then a simple induction of (3.8) using the bound (|3.7[) yields the lower 



bound 



log 4 



v{0,0)>MO,0)>a(^^\"^' minv^:f{s), s€(0,r]. 
Vlog^y K 

A straightforward computation verifies that /(s) — >■ and f{s)/s -^ 00 as s \, 0+ since 
a < 1. n 
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3.14 Corollary (Continuous expansion of {v < 0}). If v is a vise, supersolution of (1.1) 
then the set {v < 0} eannot expand discontinuously, i.e. 



{w < 0, t < s} == {w < 0, t< s}, 
for any s. 
Proof. Suppose that the claim is not true, i.e. there is a point 



(?, r) e {w < 0, t<T}\{v<0, t < t} 
That means that there is p > such that 

Sp(c,T-p)n{w<o, t<T} = 

and thus v > in Sp(^, r — p). Lemma 



3.13 



then yields w(^,t) > 0, a contradiction. 



a 



3.15 Remark. Note that the set {u > 0} of a subsolution u can expand discontinuously when a 
part of the set {u < 0} is pinched off by a collision of two "fingers" of {u > 0}. 

3.16 Lemma. Let u > be a bounded subsolution of the parabolic problem in a parabolic 
neighborhood of S, with S = Sr(C, f) for some (^, a) € M" x M and r € (0, 1). Assume that 
u = Q on E. Then there exists e > and g G C^([0,r + e)) with g — g' — on [0,r] such 
that 

< u{x,t) < g{\x- ^1) for |a; - ^| < r + e, t e [a, a + r]. 

r]. For given £ > 0, t; e (0, e), we define the 



Proof. Let us choose (y, s) G Br{x) x [a, a 
barrier 



^Criix^t) 



AM 



{4nAt 



■v), 



E, 



(SeV2 X (-^' 0]) n {^e,n{x, t) > 0}. 



Note that: 



(i) Due to (1.6), f/'g ^ is a strict supersolution of the parabolic problem in E^ ^ as long as 
^'/' < 31^' and ^e.v > 2M on dB,,,. x [-^^,0]. 



(ii) Ve,j? > if and only if i > -j^ or |a;| > ^/-AnAt - rj. 
(iii) ^/^^ > V-47iAt when te{-^,0) and -^^^^ > e^/^ for small e > 0. 
(i)-(iii) verify that for small e > and all 77 £ (0, e), we have the ordering 
u{x,t) < ip^^^{x-y,t~ s) on dpE^.,, + {y,s) 



and hence by the parabolic comparison (Prop. A. 2) the ordering holds in E^_ri 
Therefore there exists a constant e > such that 

u{x, t) < gi\x - ^1) := inf ■(/'e.r,(a; - C, 0) 

r;>0 
for (x, t) e Br+e X [cr, CT + r]. 



(y>s) 



a 
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3.3 Z and W cross 



We proceed with the the proof of the comparison principle (Theorem 3.11 for the regularized 



solutions Z and VF, in place of u and v. To argue by contradiction, we investigate the 
situation when Z and W cross in Q^, i-e. there is a finite first crossing time ioj defined by 



io :=sup{T I Z(-,t) < W{-,t) for < t < r}. 



(3.9) 



Now since {W^ < 0} cannot expand discontinuously due to Lemma 3.14 



we can prove 



that a certain ordering of level sets of W and Z is preserved up to the crossing time t^. 



3.17 Lemma. Ltt Z and W be the regularized solutions defined in (3.2) and tg be the 
crossing time defined in (|3.9[). Then 



In particular, 



{z > o},„ n {w^ < 0},^ ^d{{z> o}J nd{{w< o}J 



mt{{Z>0}Jnmt{{W<0}J=t 



(3.10) 



(3.11) 



Proof. First observe that {W < 0}^ C {Z < 0}^ for t < to (or equivalently {Z > 0}^ C 
{W > 0}, for t < to). 

Step 1. We claim that int {W < 0}^^ C {Z < 0}^^ 



Pick any ^ G int {W < 0}^ . Due to continuous expansion of {W < 0} (Corollary 3.141, 
there exists S > such that Bs{^) C int {W < 0}^ for a short time before to, t £ [to — d, to\ . 
Lemma 3.12 yields that Bs{£,) is connected to dflr with non-positive values of W (and thus 
negative Z) through a wide "trunk" 

G = (7([0,l]) + S,/2)nO,, 

and Bs{£,) C G C {W < 0}^ . The continuous expansion of {W < 0} again guarantees that 
G can be chosen so that G C {W < 0}^ for a short time before ig- Finally, we recall that 
Z < W in G X [0, to). Therefore we can construct a strict classical supersolution of ( 1.1 ) up 
to the time to that will stay above Z and is negative in Bs{£^). Indeed, we solve the elliptic 
problem in G with zero on dG D fir and negative data on dG D dVir. Since G has a uniform 
interior ball property at points dG n fir, we can proceed as in the proof of Lemma [4?T| We 
conclude Z(^, to) < 0. 

Step 2. Since {Z > 0} n Qr = {Z < 0}" n Q^, we have 

int ({W^ < 0},J n {Z > 0},„ = 0, (3.12) 

which together with Lemma |3.11| also gives 

int {{Z > 0}J n{W< 0},^ = 0. (3.13) 



Step 3. Now (3.12) clearly implies 



{Z > o},„ n {M^ < 0},„ C{W< 0},„ \ int {{W < 0}J ^d{{w< 0}J 



Using (3.13) symmetrically concludes the proof of the lemma. 



a 
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3.18 Corollary. Let Z, W and to be defined as above and let 

ee{z>0},,n{M^<0},^. 

Then there is a unique unit vector v such that v is the unit outer normal to {Z > O}^ at ^ 
and the unit outer normal vector to {W > 0}^ at ^. 



Proof. Due to Lemma 3.17 ^ e d {{Z > 0}^J n d {{W < Oj^J. Moreover there are dual 
points 

i(u,<Ju) e n"(P) C d{u > 0} and (e.,a,) G n„(P) C d{v < 0} 



(see Proposition|3.10l 



Let Bu = int I ^r(C«j '''«)*_, ) ^-nd B^ = int I Sr(^t,,crt,) ). Due to the definition of 

Er , Bu and By are balls in K" , centered at ^„ and ^„ , respectively, of radius greater than 
or equal to r. Observe that i?„ C int {W > 0}^ and B^ C int {Z < 0}^ . Therefore 

Bu n By C (int {W > OJtJ n (int {Z < 0}J = 0, 

while Bu n By = {£_}. This is true for an arbitrary choice of dual points. Therefore the 
outer unit normal of v of i?„ (resp. inner unit normal of By) at f is uniquely determined, 
and can be taken as 



le-Cni i^.-er 



a 



3.4 Finite speed of expansion 

Our goal in this section is to use the ordering of the support to prove the ordering of the 
functions Z and W at the contact time i = to- This needs a careful analysis since Z and W 
are merely semi-continuous. 



From Lemma 3.3 we know that Z_|_ is a parabolic subsolution, and —W^ is an elliptic 
supersolution for each time, and also W^ is a parabolic supersolution in {W > 0} (open set) 
and Z is an elliptic subsolution for each time in {Z < 0} (open set). Therefore we can invoke 



the standard comparison principle (Propositions |A.l" and A. 2 1 in the open sets {W > 0} and 



{Z < 0}(, if we know that the functions are ordered on the boundaries of these sets. This is 
not completely obvious, however, and we have to pay special attention to the situation at 
the contact point. 

3.19 Lemma (Finite speed of expansion at the contact point). Let Z , W be the regulariza- 



tions and to be the crossing time as defined in (3.9). Then for any ^ € {Z > 0}^ n{VF < 0}^ 



n"(e,to)un„(e,io)c9iS,(e,to). 

Proof. We split the proof into a number of shorter steps. Here we use the simplified notation 
S = Sr(^,to) (and analogously for H" and n„). 



Step 1. n„n 9^2 = 0. 

For the sake of the argument, suppose that this does not hold and there indeed is a 
point P' = {C, a') e n,„ n 9^S. Note that a' = to + r. We observe that w > in S due to 
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Figure 4: Construction of a test function tp in Lemma 3.19 



(3.10) and Proposition 3.10 |iii), and thus we can apply Lemma 



consequence is |^' — S,\ = r, i.e. 

P' (^dDr{S,M + r 

Let V be the unit normal from Corollary 



3.13 



for V on :i. The first 



in Prop. B.3) and set 



3.18 



We choose po = min(|,/5c) {pc is defined 



C, = ^' — poiy, a = a', a = 2, b 



-1 and a) = 0. 



Proposition B.3 provides us with a radially symmetric subsolution of ( |1.1[ ) with parameters 
d, b and w. We denote tp its translation by (C,^) which is defined on the cylinder 



K = {po - £ < \x - C\ < Po + e, te[&-e, a]}, 



(3.14) 



for some e > 0. The situation is depicted in Figure |4J which is a part of the larger picture 
in Figure [SJ 

Let us define ip = fj,ip, where /i > is picked large enough to ensure that min;j^(^ < 
min^j^t;. Furthermore, since v{x,t) > f{r— \x — ^\) for |a; — ^| < r, with f{s)/s — > cx) as 
s — > 0+, there exists ei > for which u > 1^9 on 

{po - £1 < |a; - CI < Po, t e[a- e, a]}. 

We shall compare v and 93 on a cylinder S, 

S := {po ~ei <\x-(\ < po + e, t e {& - e, &]} c K. 

So far we have shown that v > ip on all of dpT, except on 

{{x, a - e) : Po < \x - (\ < Po + e} C E, 

where clearly p < while v > 0. 

Therefore p < v on dpYi while p > v at P' £ Y,, and we arrive at a contradiction since 



(^ is a strict classical subsolution of (1.1) on S 



Step 2. n" n S^S 
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The proof of step 2 is similar to that of step 1. Choose P' E H" n 9^^. Due to (3.11 1, 
Lemma 3.12 and Hopf's lemma (Prop.[AT3), we have 



li^infZ^il±i^=2c>0. 

h\0+ h 

Therefore there is ft-o > for which 

u{x, t) < Z{^ + vh, to) < -ch < ior h<ho, {x, t) G ^^(C + yh, to). (3.15) 

Proceeding as in the proof of step 1, we have (p, e € (0, /iq) and K, this time with C, — C+Poi^, 
a = a' . In contrast with step 1, we define (p = —fJ-'P, and (3.15) lets us choose /i > small 
for which ip > u on 

{\x -C\= po-e, t e [a - e,d-]} n {po - e < \x - C\ < Pa, t = a- e}. 

Now we find £i > such that ip > u on {\x — CI — Po + ^ii ^ G ['5' — £, o']}, this is again 
possible thanks to Lemma [3.16| applied to w on S. 
Finally, define 

E := {po - e < \x - C\ < Po + £i, t e {a - e, a]}. 

Since clearly (/5 > and u < on {po < \x ~ (^\ < po + £, t — a — s}, we succeeded in 
showing that p > u on dpT, while Q = p < u a,t P' ^ T,^ a, contradiction. 



Step 3. 11" n 5j_S ^ then l it, H d '^E = (and also when if we swap 11" and ![,„) 

Indeed, this follows from (3.11) and the definition of to in (3.9), which together yield 



that for any P„ € 11" and P^ G 11^ we have 

s,(PJ n s,(p,) c^{t< to} = (e, to) = P. 

step 4. To finish the proof, we simply realize that steps 2 and 3 used together imply that 
n" n Sj^S = 0, and similarly steps 1 and 3 imply n„ n 9j^S = 0. D 



3.20 Remark. Lemma [3?T9| shows that the situation at any P = (5, to) '^ {Z > 0}r\{W < 0} looks 
like Figure |5] 



Now we have enough regularity to show the following: 

3.21 Lemma. For Z . W , U) defined above, we have 

Z^W^^O on {Z>0}j^n{IF<0}t^. 

Proof. We will only show this result for Z. A similar, simpler argument applies to W as 
well. See also |KP[ Lemma 3.6]. 

Let M = 2 ma.X[t<to} Z < oo and P = (^, to) e {Z > 0} n {W < 0} We can also choose 
Py e n„(P). Let v be the unit normal from Corollary 3.18 Lemma [3.19 guarantees that 
there exists rriy € M such that {iy,my) is an interior normal of d'E{Py) at P, see Figure [5] 
This is why we call Lemma |3.19| "finite speed of expansion" . 

We recall that Z < W < in S(P„) n {t < to}. Furthermore, the argument from step 



Z < in 5(P„) n{t = to} as well. 



2 of the proof of Lemma 3.19 using (3.11 ), Lemma 3.12 and Proposition A. 3 verifies that 
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(e,to)^^^ 



\{M^<0} 
Figure 5: Situation at P = (^, to) e {Z > 0} n {W < 0}. 
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X • I^ 



ft, Pq — LOT 

Figure 6: Construction of a test function ip in Lemma |3.21 



to - r 



As in the proof of Lemma 3.19 we can construct a strict classical supersolution of 



(1.1) with the help of Proposition B.3 Let us again choose po small enough so that for 



a = 2, b = —1, Lu — max{0, 2mi,}, C = C + Po^j there is a strict classical supersolution 
—(p{- — (,■ — to), defined on the set 

K ^ {po - e < \x - C\ < Po + e, te{to-e, to]}, 

for some e > 0. 

For T > and h = t^ , let us define 

Sr = {po <\x-C,\-{t- to)Ld < po + h, t e (to - T, to]}, 

see Figure [6j 

We can take r small enough so that 

(a) Sr C K, 

(b) S, n {t = to - r} C S(P')- 
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Now we find /i > large enough so that 

-fiip>2M on {\x\= po + h + tu), t e [~T,0]}. 



For 77 > set (pri{x, t) — —iitp[x — (^ — rjv, t~to). Since Z is a subsolution of (1.1 ) and Z < ip 
on dpY. + [rjv, 0) for small 77, we also have Z < (ys in E + {rjv, 0). We conclude by sending 



Z{£„t(i) < lim iprji£„to) = ^oi^,to) = 0. 



a 



Lemmas 3.17 and 3.21 have the following consequence: 
3.22 Corollary. Let Z, W and to be as above. Then 

Z <W att^to. 
Furthermore, 

{z =. 0},^, n{w^ o},„ ^ 0, 

Proof. We apply the comparison principle for elliptic equations (Prop. |A.1[ ) and compare 
Z(-,io) and W{-,tQ) on the open set {Z < 0}^ . Indeed, we can set Z = on d{Z < 0} to 

ensure that Z < Q. The modified Z{-,to) is clearly in USC{{Z < 0}^ ) and a subsolution 

of the eUiptic problem in int {Z < 0}^ (= {Z < 0}^ by Lemma |3. 11 1. Then Z < W on 



d [{Z < 0}j ) due to Lemma 3.21 and the elliptic comparison applies 



The same can be done using the comparison principle for the parabolic equation (Propo- 



sition [AT2| ) in the parabolic neighborhood {W > 0, t < to}. 

Finally, if {Z = 0}^^ n {W = 0}^^^ = 0, then in the view of Lemma 
{Z > 0, t<to}n{W <0, t< to} = 0. 



3.17 



and Lemma 



3.21 



Since these sets are compact, they have positive distance and the comparison in the first 
part of the proof yields Z < W at t = to, a. contradiction with the definition of to and 

z-w e use. a 



3.5 Ordering of gradients at the contact point 

Here we will show that, at a contact point P, the "gradients" of Z (resp. W) follow the fiux 
ordering given in Definition |2.2[ iv) (resp. its supersolution counterpart). As mentioned in 
the beginning of section [3j such ordering would readily yield a contradiction with the fact 
that Z crosses W from below at P. 

3.23 Lemma. Let Z , W and to be as defined above, let 

P = (^,to) G {Z > 0} n {W < 0} be a contact point at time to, and let v be the unique 

spatial unit normal vector at P obtained in Corollary\3.18\ Then 



lim inf 

h\0+ 



Z{^- hv,to) 



lim sup 

h\0+ 



wii- 



h 
hv,to) 



liminf^(i±^>0. 



lim sup 



W{i + hv,to) 



<0. 
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Proof. We prove the result for Z. The proof for W is similar. 

Let us again write S — Sr(P) = Er{£,,to) throughout the proof. Denote 



lim inf 

h\0+ 



Z{^ - hv, io) 



h :— lim inf 



Z{£, + hv,ta) 



Note that 6 < < a. A simple barrier argument, following the one in the proof of 
Lemma |3 . 2 1 1 and taking advantage of the room provided by Lemma |3.19[ shows that a < oo. 
Furthermore, the inequality a + 6 > is satisfied trivially if 6 = 0. Thus we may replace h 
by —2a if necessary, and assume that — oo < & < 0. 

Suppose that the result does not hold, in other words, a + 6 < 0. In that case, we set 



K := s^ > 0, and we note that 



b< 



a < K, 
Let us choose r/ > such that 

a + 3ry < K, 6 + 3?] < —k. 

We shall use this to construct a barrier that crosses u from above at an arbitrary fixed 



since 



point Pu = iC , cr') G n"(P), yielding a contradiction. As in the proof of Lemma 3.21 

Pu € d^S thanks to Lemma 3.19 there is ?ti„ e M such that {v, rriy) G M" x M is an interior 

normal vector to 9S at Pu (see Figure Isl) . 



Following the proof of Lemma 3.19 we can construct a radial supersolution of (1.1) 



with the help of Proposition |B.3[ Let us again choose po small so that for a = a + 3?? and 
h = b + Srj, uj = max(0,2m), C = C' + Po^^i and ^ = cr', there is ip — ~(p, a strict classical 
supersolution of (1.1) translated by (Ci^), defined on K as in (3.14), and ip ~ a.t Pu- 



From the definition of a and b, and Z, for every r e (0, 1) there are hi,h2 € (0, r^) such 



that 



sup u ^ Z{^ — hiv,to) < {a + 'r])hi, 



sup u 



Z{i + h2V,to) <{b + ri)h2 



(3.16) 
(3.17) 



For a given choice of r, /ii, /12, let us define 

Sr := {/Oo -h2 <\x- C,\- [t~ a)u) < po + hi, t e (a - t, a]}, 

see Figure [4] for a sketch of a similar construction. We shall choose r £ (0, 1) small enough 
so that the following holds: 

(a) S, C K, 

(b) p > {a + r])hi on A := {\x — C\ — {t - ajib = po + hi, t e [a ~ r, a]}, 

(c) p> {b + T])h2 on dpJ^r \ A, 

(d) AcEr{C-hii^,tQ), 

(e) dpJ:r\AcEri^ + h2V,to). 
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Indeed, (p satisfies (b) and (c) for small t due to its smoothness in the positive and negative 
phases, and (a), (d) and (e) are a consequence of the choice of Cj and the definition of E^-. 
Since ip > u on dpUr (thanks to (b) + (d)+( pl6l ) and (c)+(e) + ( |3l7| ) while = (^ = m at 
P„ e St, we get a contradiction. D 



3.6 Proof of Theorem l3J] 

Now we are ready to prove Theorem |3.1| The proof proceeds by showing the comparison 



for the regularizations Z and W defined in (|3.2|. We choose r > such that Z < W on 
dpQr (Proposition 3 



Suppose that there is a point in Q where u > v. Since u < Z and W < v, we see that 



there must be a point in Qr where Z > W. Then to in (3.9) is finite. We recall that Z <W 
at t = to by Corollary |3.22 



Corollary |3.22| also guarantees the existence of a contact point at t = io : 
P = (C, to) e d{Z > 0} n {T4^ < 0} with Z{P) = W{P) = 0. 
At the point P, the boundaries d{Z > 0} and d{W < 0} are C^'^, in the sense that they 



have space-time and space balls from both sides (Proposition 3.10 iii))- Let v be the unique 
unit normal vector to the space balls at ^ given by Corollary |3.18| 

Set n' to be the connected component of the open set {Z < 0}^ that contains £_ on its 
boundary. Lemma |3. 12| guarantees that 

Z{;to)^W{-,to)onn'. 

We recall that Z{-,to) is a subsolution of the elliptic problem on H.' and —W~{-,tQ) is a 
supersolution of the elliptic problem (Lemma 3.3). We can apply the elliptic Hopf's lemma 
(Proposition A. 3) to Z{-,to) and W{-,to) on il' at ^, which yields 



hminf Z^i^ll}^ > H^,,f-^(e + ^Mo) 

h\Q+ h h\Q+ h 



Finally, the "weak gradients" are ordered at ^ by Lemma [3. 23| which leads to a contra- 
diction (we write Z{-) instead of Z{-,to), all hmits are h \ 0-f): 

,. .,^(C-H^r -^(e + H^r . .-^(e + H 

hm mi ; > lim sup > lim mi ; 

h h h 

Hopf .^-VK(e + H^,. w^(C-M 

> hm mt > lim sup , 



> liminf 



h 

Z{i - hv) 



Therefore to cannot be finite and we conclude that u < v in Q. This finishes the proof of 
Theorem O D 
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3.7 



Remarks on the proof of Theorem 3.1 for divergence-form op- 
erators 



The arguments in the proof of Theorem 3.1 also hold for the divergence- form operator F 



given by (1.4). The only difference lies in the references on the properties of solutions of 



the parabolic and elliptic problems, which have been used throughout section [3] and which 
are proved or referred to in Appendix [Xj as well as the specific barriers constructed in 
Appendix |B.1[ 



Since the barriers for the divergence-form operator are constructed in Appendix |B.2[ here 
we only point out the references for the regularity properties of solutions of the parabolic 



and elliptic problems (in the sense of Definition 3.2). More precisely, we used the parabolic 



Harnack inequality (Proposition |A.4[ ) and the elliptic Hopf's lemma (Proposition A. 3). For 
a divergence-form operator F given by (1.4), Proposition A. 4 is shown in [LSU[ Chapter V] 
and Proposition A. 3 is shown in IDDS . 



Having all of the above properties and test functions, one can proceed as in section [3] to 
prove the following: 



3.24 Theorem. Theorem 3.1 holds for F given in (1.4) 



4 Existence and stability of viscosity solutions 

Let F and 17 be given as before. In this section we show the existence of viscosity solutions 
of ( 1.1 1 via the Perron's method. We say that a continuous function uq : i^ ^- K is in T' if 

F{D'^uo, Duq, uo) = in {uq < 0}. 



the following conditions hold: 

uq = -1 on dil and - F{D'^uq, Dup,, uq) = in {uq < 0}. (4.1) 

r(wo) := d{uo > 0} = d{uo < 0} is locally a C^'^-graph. (4.2) 

Note that (14. 21) and the boundedness of il guarantee that there exists Rq > such that 



d{uo > 0} has both an interior and an exterior ball of radius Rq at each point. 

First we ensure, by constructing suitable barriers, that the solutions with initial data 
Uq €z V evolve continuously at the initial time. 

4.1 Lemma. Suppose uq £ V. Then there exist at least one viscosity subsolution U and 



one supersolution V of (1.1) with initial data uq. Moreover, there exists a constant C > 



and a small time to > depending only on the maximum of Uq, N , d and Rq such that for 
< t < ip we have 

d(x, F), d(y, F) < t^l^ for x e d{U{-, t) > 0} and y £ d{V{-, t) < 0}. 

Proof. 1. First we construct a supersolution V. Let us set 

0+{t) := {x : d{x,{uQ > 0}) < i^/^}, 

and let V{x,t) solve 

Vt-F{D'^V,DV,V) >0 

V 







-F{D^V, DV, V)>0 



0+{t) 

r{t) -.^do+it) 
o-{t) ■.= n\o+{t) 
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with the initial and lateral boundary data uq and —1, respectively. Note that, due to the 



regularity assumption (4.2) on F, the positive phase C+(i) has the exterior ball property 
with radius Rq — t^'^ for small time t. In particular, it follows from the uniform ellipticity 
of the operator in the negative phase that there exists to > such that we have \DV\ > cq 
on r(i) for < i < to, where cq > is independent of t. 

Now a barrier argument, which takes advantage of the barrier constructed in the proof 



of Lemma 3.16, yields the existence of a function c{t) such that 

\DV+\ < c{t) with c{t) ^- as i ^ oo. 
Hence we can choose Iq sufficiently small such that |-DT^^| < |£^V^~| on T{t) for < 



t < to, and therefore T^ is a viscosity supersolution of (l.ll for < t < to- For t > 



to we take V{-,t) — ijj, where ipix) is a viscosity supersolution of the elliptic problem 
—F{D'^iIj, Dip, ili) = in rj with the boundary data —1, and -0 > V{-, s) for < s < io- 

2. The construction of a subsolution U is similar: we replace C+(i) in the construction 
of y by 

0+{t) := {x : d{x,{uo < 0}) > t^^^}. 

It follows from a barrier argument in the proof of Lemma |3.13| that 

\DU+\ > C{t) on F(f), with C{t) -> oo as i ^ 0. 

On the other hand due to the exterior ball property l-DC/^j < Cq on F(t), where Co is 
independent of time. Hence we can choose to sufficiently small such that |-Df7^| > |-DC/~| 
on F{t), and thus U is a viscosity subsolution of ( |1.1[ ) for < t < to- For t > to we take 
C+(t):=0. a 

Using V and U constructed above, the minimal solution u can be constructed as 



u — inf{z : viscosity supersolution of (1.1 ) such that U < z}, 

and the maximal solution u can be constructed as 

u — snplw : viscosity subsolution of ( |1.1[ ) such that w < V}. 

4.2 Theorem, u and u are viscosity solutions of (|1.1[ ) with initial data Uo, and u = u* , 
u = u^. Moreover, if v is a viscosity solution of (1.1) with initial data Uq then u <v < u. 



Proof. The proof follows from standard arguments in Perron's method (see for example the 
proof of Theorem 1.2 in K3 . Also see GIL section 4]). D 



Lastly we show stability properties of the maximal (and minimal) solutions under per- 
turbation of initial data. 

For a family of functions {u'^}e>Oi let us define the half-relaxed limits 



and 



limsup* u'^{x,t) 

£->-0 



liminf* u'^ (x, t) 



limsup u'^{y,s) 

e-i-O 



liminf u'^(y,s). 

{y,s)^{x,t) 
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4.3 Theorem. Let u he the maximal viscosity solution of ( |1.1[ ) with initial data uq. Let ' 
he a sequence of Junctions decreasing in e ~> such that 

(i) uq < Uq and Uq converges uniformly to uq as e — >■ 0, 



(ii) Uq and F^ :— d{uQ > 0} = d{uQ < 0} satisfy Assumptions (4.1) and (4.2) uniformly 
in e > 0, 

(Hi) r^ converges uniformly to r(uo) with respect to the Hausdorff distance. 

Let u^ be a viscosity solution with initial data Uq. Then m^ converges to u in the following 
sense: both Ui := limsup*g_j.o u^ and U2 '■— liminf*£_>.o m"^ satisfy 

(Ui)* = u and {Ui)^ = u*. 

Corresponding results hold for the minimal solutions, if one replaces the condition uq < Uq 
in (i) with Uq < uq, and Uq is increasing as e — >■ 0. 

4.4 Remark. Note that the convergence in the sense of the above theorem is optimal for semi- 
continuous functions. 



Proof. 1. First observe that Theorem |3 . 1| yields that 

u* < [u^)^ for any £ > 0. (4.3) 

On the other hand, due to the standard stability property of viscosity solutions, 



limsup ^^QU^{y, s) is a subsolution of (1.1). Hence by a barrier argument (one may use 



the supersolution V^ constructed as V in the proof of Lemma 4.1 but for the initial data 



Uq) we can show that 



limsup* u^<K (4.4) 



Therefore by definition of u we have 



limsup* u'^ < u. (4.5) 

£->-0 



Putting (4.3) and (4.5) together, we conclude. D 



From the above theorem we have the following "regularity" information on the minimal 
and maximal viscosity solutions. 

4.5 Corollary. 

(u*)* — u and (u*)* — u. (4-6) 

5 Uniqueness properties 

Recall that our comparison principle. Theorem 1 3. 1[ requires strictly separated initial data. 
We saw in the previous section that the theorem yields existence (and uniqueness, by defini- 
tion) of maximal and minimal viscosity solutions. In this section we will discuss uniqueness 
properties of general viscosity solutions. 
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5.1 Viscosity solutions coincide with regular weak solutions 

In this section we consider a linear uniformly elliptic operator F of the form F{M,p, z) = 
tr M. We show that in this case the viscosity solutions of ( |1.1[ ) are precisely the regular weak 
solutions introduced in AL . As a corollary it follows that viscosity solutions are unique for 
this class of operators. 

For completeness of the paper we revisit the notion of regular weak solutions and its 



properties. We refer to Section 1.4 and Definition 2.1 of AL for the definition of weak solu- 
tions, subsolutions and supersolutions of (1.1 ), which are as usual defined in LF'{Q, T; H^{n)) 
via integration by parts. 



5.1 Theorem ( [AL| , Theorem 2.2: simplified version). Suppose that 

F{M,p,z) = trM, 

and let Q = Q X (0,T], where CI is a bounded domain with Lipschitz boundary. Suppose that 
ui is a weak subsolution and U2 a weak supersolution of (1.1), with initial data b^, fojj o.iT'd 
boundary data uf , u^, respectively, such that 

dtb{ui), dtb{u2) e L\Q). (5.1) 

If b1 < 62 a.e. in J7 and u^ < u^ a.e. on OlQ then ui < U2 «-e. in Q. 

We call u a regular weak solution if u is a weak solution and satisfies the additional 



regularity (5.1). It is shown in |AL| (Theorem 2.3) that for linear F there exists a regular 
weak solution with initial data uq G V. 

Observing that the test functions U and V constructed in the proof of Lemma |4.1| are 



respectively a classical subsolution and a classical supersolution of (1.1) for < t < to, 
Theorem |5 . 1| yields the following: 



5.2 Lemma. Let uq G H^{Q,) satisfy assumptions (4.1)-(4.2|. Then any regular weak 
solution u of (1.1 1 with initial data uq uniformly converges to uq as f — )• 0. 



Furthermore, the following holds due to the estimates derived in the proof of Theorem 



2.3. in AL 



5.3 Lemma. Regular weak solutions are stable under perturbation of initial data in H^(fl). 
More precisely if the sequence of initial data u'^ converges to Uq in H^(ri), then the corre- 
sponding regular weak solutions Un with respect to the initial data m" converges to the regular 
weak solution of u with initial data uq in L^{0,T; H^{fl)). 

Given a weak solution u € L°^[Q), we want to find a representative v = u a.e. that 
is a suitable candidate for a viscosity solution. We shall find it using a weaker notion of 
semi-continuous envelopes. To this end, let us state an elementary lemma. 

5.4 Lemma. Let E C M" such that \E\^ B&{x)\ > for all x e E and S > 0. If f < g a.e. 
in E, then 



f> 



*.E 



< 



9 
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Given a function u G L°°{Q), where Q = fi x (0, T], define tlie essential semi- continuous 
envelopes 

u^{x,t) = inf v*''^{x,t) = inf ess sup u, 

v=u a.c. r>0 „ , ,^„7^ 

U(^(x,t) — sup v^ Q{x,t) — sup ess inf_u. 

«=u a.o. r>0 -Br(x,t)nQ 



Here v*''^ and u^ g are as defined in ( |f.9[ ). 

It follows from the definition that u^ e USC{Q) and ti^ G LSC{Q), and due to Lemma 
15.41 we have 

We claim that u^ < u < u^ a.e. To see this, we have for almost every P = (a;, t) by the 
Lebesgue's differentiation theorem 

u^(P) = inf v*{P) ~ inf inf sup v 

v—u a.c. iJ— It a.c. r>0 d ^ p-j 

= inf inf sup v > lim — -- / u = u{P), a.e. P. 

r>0 v=u a..c. g^f^pj r^Q \Br(P)\ J B,.(P) 

Next we introduce a candidate function for the viscosity solution. 

5.5 Lemma. Given a function u € L°°{Q), define 

V = max(min(w, w*), Uo). 
Then v = u a.e., v■^, — u^ and v* = u^ . 

Proof. Observe that u^ < w < u* since u^ < u^ . In fact, since u^ < m < u^ a.e we also have 
that V = u a.e. Now u^ £ USC implies v* < u^ and the definition of u^ as the infimum 
yields u^ < v* . Similarly U(^ — v^. D 

Now we are ready to state the main result in this subsection: 

5.6 Theorem. Let V be the class of regular initial data as defined in Section p] and let 



be the unique regular weak solution with initial data uq £ V . Then v as given in Lemma 5.5 



is a viscosity solution of (l.I) for the initial data uq 



Proof. The proof is based on the local comparison principle that regular weak solutions sat- 



isfy (Theorem 5.1 1. First note that v uniformly converges to uq at t = due to Theorem 5.1 



and the barrier arguments using the barriers V and U constructed in Lemma |4.1 



We only show that w* is a viscosity supersolution; parallel arguments apply to the subso- 
lution part. Let </) be a strict classical subsolution in a cylindrical domain Q' := fl' x (ii, t2] C 
Q such that < u* on dpQ' . 

Our goal is to show that cj) < v^, in Q'. Thus assume that this does not hold and there 
exists (a;o,io) G Q' with > v* at (xojio)- 

We claim that we can perturb (j) into a strict classical subsolution (f> with > at (a;o, ^o) 
and < u* on dpQ' . To this end, let 6 € C^(M"), 6' > 0, be the standard smooth moUifier 
with support B^. For e,r] > let 6^^^{x) := e9{i]^^x) and define <p as 
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• if 0(a;o,to) = 0: tf>{x,t) = (j){x + 9g_,j{x)h',t), where v is the unit outer normal of 
{(/) > 0}t^ at a;o, 

• if (l>{xo,to) 7^ 0: 4>{x, t) = (j){x, t) + d^^,^{x - xq)- 

If £ and rj are chosen small enough, a straightforward differentiation shows that the pertur- 
bation (j) has the required properties. 

There is also S,r] > such that (p + S < v^, on an 77-neighborhood of dpQ' , i.e. on the set 

{(a;,i) eQ' ■.dist{{x,t),dpQ') < 77}. 

Finally, v^ < u a.e. in Q' by definition. We conclude that: 

(i) (/>{■, t) < u{-,t) on dfl' for a.e. t e (^1,^2] in the sense of trace on H^{n'); 



(ii) due to regularity ([sl]), 6(u) G C{ti,t2;L^{n')) (see [e) §5.9.2 Theorem 3] or |GGZ 
Chapter IV, Theorem 1.17]) and therefore there is a unique 6q G L^(51') such that 
b{u{t)) — 7- 6q strongly in i^(^') as i — >■ ti, and (f){-,ti) < b'g a.e. on ^' . 

A simple computation shows that (p is a, regular weak subsolution in the sense of 



Definition 2.1]. Therefore, by Theorem 5.1 we have 4> < u a.e. in Q' and hence <f> < 



AL 



< Wo = w* at (a;o,to), a contradiction. Therefore w* is a viscosity supersolution of (1.1) 



due to Remark 12.51 D 

Theorem |4.3| and Lemma [5 . 3 1 yield the following: 

5.7 Corollary. Let F{M,p,z) — tr M and let uq £ V. Then there exists a unique viscosity 
solution V of ( |1.1[ ) with initial data Uq, which coincides almost everywhere in Q = V,x {0, c»] 
with the regular weak solution with initial data uq. 

5.2 Discussion of uniqueness for nonlinear operators 

At the moment, the uniqueness results for fully nonlinear operators are unfortunately quite 
limited to initial data and domains under special geometries. Here we illustrate an example. 
The authors suspect that uniqueness holds in most settings, but it remains an open question. 
When 

b{u)^u+, F ^ F{M) with F{fiM)^ fiF{M) for ^i>0, (5.2) 

the following holds. 

Let Fq := d{uo > 0} = d{uQ < 0} be a locally C^'^ graph, and suppose the initial data 
Uq satisfies one of the following: 



(A) (Initially decreasing: only possible) uq is C^ in {uq > 0} with 

-F{D^uo) > in {uq > 0} and \Du+\ < \Du^\ on Fq; 

(B) (Star-shaped) g — —1, 57 is a star-shaped domain, and uo((l + e)x) < uo{x) for all 
X G {1 + e)^^fl except at x = for any e > 0. 



5.8 Theorem. Suppose that (A) or (B) holds. Then for the operator of the form (5.2) 



there exists a unique viscosity solution of (1.1) with initial data uq. 
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Proof. 1. In the case of (A), let us define ^0(2;) := {uq — e)+ in {uq > e} and Uq be the 
solution of 

-F{D^u) = in {uo < e} 

with the boundary data g. Since uq is C^ and D{uq) changes continuously with respect to 
the change of the domain (due to the ellipticity of F and the regularity of the domain) , if e 
is sufficiently small we have |D(wo)'''| < \D{uo)^\. 

Now let us define V{x,t) := Uq{x). Then for small time, V{x,t) is a viscosity supersolu- 



tion of (I.l ). It follows then that, for sufficiently small e > 0, 

u{x,e) < u(a;,0) in i7, (5-3) 

and there exists a constant q — >■ 0"*" as e — >■ such that (1 + Cg)u(a;,e) < u{x, 0) in fi. Let 



u and V be two viscosity solutions of (1.1) with initial data uq. In case of (A) we perturb 



by (1 + cos)u{x, t + e) for given e > and apply (5.3 1 as well as Theorem 3.1 to show that 
u < V. 

2. In case of (B), first note that from the star-shapedness assumption it follows that 
wo(0) > 0. Hence for given e > 

(1 + c,)u{{l + e)x, (1 + eft) < u in (1 + e^^n x [0, 00). 

Since e is arbitrary, we conclude that u < v. D 

6 Approximation by uniformly parabolic problems 

In this section we show that our problem ( |1.1[ ) can be obtained as a singular limit of uniformly 
parabolic problems. Consider bn € C^(M) to be a smooth approximation of b(s) = s-|_ with 
the following properties: 

(a) < fe;, < 1 on M, 

(b) bn —T' b locally uniformly on M, 

(c) b[^ — > locally uniformly on (—00, 0], 

(d) b'n — > 1 locally uniformly on (0, 00). 
For example, one can choose 

1 e" + e"^'* 
bn(s) — -^ log . 

6.1 Remark. The interval in condition (c) includes 0. This moves the irregularity of b^ into the 
positive phase and thus simplifies the proof of Proposition |6.2| but it is not an essential assumption. 

For n > 1, let Un be the unique solution of the smooth parabolic problem 



I K{un)t - F{D'^Un, DUn, Un) = 0, in Q := ft X (0, oo), 
[un^g on dpQ. 

Here we assume that <? is a bounded Lipschitz continuous function on dpQ. 
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(6.1) 



Due to the maximum principle, m„'s are bounded uniformly in n and therefore problem 
(|6.1) is uniformly parabolic for each n. Note that m„ is C^'" in Q for each n due to LSU 



(for F given as in (1.4)) and W2| (for nonlinear F). Also note that one cannot directly 



use the results from [GIL] since the operator does not converge uniformly. Indeed, b' is 
discontinuous in our setting. 
Define the half-relaxed limits 

uj := limsup*M„, w := liminf*u„. 

n — ^cx2 ^ ^^^ 

6.2 Proposition, w and lo are respectively a viscosity subsolution and a viscosity superso- 
lution of (|1.1[) in Q. 



Proof. Let us denote u = UJ. The goal is to show that w is a viscosity subsolution of (1.1 ). 
1. Suppose not. By definition, there exists S > and a strict classical supersolution 
of ( |1.1[ ) in a parabolic neighborhood E — U Ci {t < t} C Q, satisfying 

(i) -F{D^(j), Del), (/)) > (5 if (/.< 0, 

(ii) (t)t - F{D^(l3, D(j), 0) > (5 if (/) > 0, 

(iii) I £>(/>+ 1 + 35 < I £>(/)- 1 on {0 = 0}, 

such that u < (f> on dpE while (j) < u at some point oi P ^ E. By increasing r if necessary. 



we 



we may assume that P G {i < r}. And by perturbing (j) as in the proof of Theorem 5.6 
may assume that in fact cj) < u ai P E E. Therefore for large n we have Un < 4> on dpE 
while u„ > at some point of E. 

2. For large n, due to the continuity of m„, we can find (x„, i„) G E with u„ — (/i = at 
{xn,tn) and w„ — < in {i < t„}. Compactness of -E allows us to select a subsequence, 
also denoted n, such that (a;„,t„) — s- (a;o,to) e E'. Since u < on SpE we conclude that 
(xo,to) G £'■ Moreover, observe that (j){xo,to) = 0, since otherwise (i)-(ii) and the locally 
uniform convergence of b'^ to b' on M \ {0} yield a contradiction. 

The regularity of u„ and the condition (iii) on cj) ensure that 0(x„, i„) y^ 0. Next suppose 
that along a subsequence we have (f){xn,tn) < 0. Since also 

[6„((^)t - F{D^(f>, D(j), (j))] (x„, i„) < 0, 

the uniform convergence of 6^ to 6' in (— oo,0] contradicts (i). 

3. Therefore we can assume that 0(x„,t„) > 0. Let us first observe that cj)f{xo,to) > 0: 
this can be verified by adding the two inequalities 

[(l)t - F{D^cl),D^,c^)] ix^,t^) > 6 
[-b'„{^)^t + F{D^, D(f>, 0)] (a;„, t„) > 0, 

and by using the fact that 6^ e (0, 1). In particular, this implies that {(/) > 0} is expanding 
at (xo,io)- 

Let V be the unit outer normal of d{(j) > 0}^ at xq. Due to the regularity of {</) = 0}, 

it is possible to choose pn > small enough such that Bpg{(o) D {(j) > 0}^^ = {2:0}, where 
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Co = xq + pqv. For small enough po, Proposition B.3 provides a strict classical supersolution 



■0 of (1.1), with parameters 

a ~ \D<j)'^{xQ,tQ)\ + S, b = — \D(f)~{xQ,tQ)\ + 6 and w = 0, 

in a neighborhood K of dBp^ x {0}. Since ui = 0, ip is independent of time and thus it is a 
classical strict supersolution of the elliptic problem in each phase. Let us define iprj{x, t) := 
ipix — Co + W^t)- Choose Si > and 77 > sufficiently small so that ipriix^t) > (j^ix^t) on 
dpY^St and S^^ d K + (C,io), where 

Ssi = Bs^{xq) X (io - (5i,to]. 

This is possible due to the gradient ordering condition (iii) of (p and due to the fact that 
0t+ > at (xo, to)- 

Note that V't; < </> at {xq, to). Therefore there is again a subsequence of n, and some other 
sequence (x„,t„) G S^^ where u„ — -0^ has a maximum zero at (a;„,i„) in E^^ n {t < t„}. 
As before, regularity of u„ ensures that ipixn, tn) 7^ and so we have 

This contradicts the fact that ip satisfies —F{D'^tp,Dip,tp){xn,tn) > 0. D 



6.3 Corollary. (For nonlinear F) Let uq £ V and Uq be as given in Theorem 4-3' o-nd 



choose a sequence Sk — > 0. Then for given T > there exists n^ — > cxd such that {uru-} 



solving (6.1 1 with initial data Ug'° converges to the maximal viscosity solution u of (1.1) 



with initial data Uq. More precisely, both 



Ui := limsup* «„,_ and U2 '■= liminf* u„j. 

satisfy 

(Ui)* = u and {Ui)■^, = u*. 

Parallel statements hold for the minimal viscosity solution. 

Proof. The main challenge in this proof is to ensure that uJ and a; assume the correct 
boundary and initial data. Let us denote AI := niax{|uo|j 1} as the uniform bound on \u„\ 
obtained by the maximum principle. 

1. We first construct barriers to control u„ near the lateral boundary data. Consider 
(p{-,t), which solves F{D'^(p,D(f,(p) = in fl for each t > with boundary data g{-,t). 
This gives a uniform bound from below on u„ at the lateral boundary d^Q. Next define 
ip{x,t) = ipix) where ip solves F{D'^ilj,Dtfj,^) = in S5 :— {x : d{x,dil) < 6} with 
boundary data g on 9J7 and M on {cc : d{x, dfl) = 5}. Here (5 > is chosen small enough so 
that f/; > g{-,0). Then m„ < </? in the set E^ and we obtain the uniform bound from above 
on u„ at the lateral boundary OlQ- 

2. To obtain estimates on w„ near the initial data, one can consider barriers which are 
defined for short time < t < io as follows. For given 5 > let V{x,t;S) be the strict 



supersolution of (1.1 1 for < t < ig, constructed as in the proof of Lemma 4.1 with the 
initial data Uq := uq + S, boundary data S on T{t) and 

0+{t) = {x : d{x, {uo > 0}) < <!/*} = {x : d{x, K > S}) < ^'^^l- 
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Note that then we have 

0+{t) ^{V>6} and 0-{t) 



n-o+{t)^{v <s}. 



We observe that Vi > in O^ {t) since V solves the elhptic equation with boundary data 5 
in O^ (t) and the set shrinks in time. Moreover, as mentioned in the proof of Lemma 4.1, 
due to the uniform elhpticity of the operator in the set 0~{t), we have 



\DV+\ < \DV-\ onr(t). 



■2) 



Now (p{x,t) := V{x,t;S) satisfies 



and 



(fit - F{D^ip, Dip, if) > -F{D^ip, Dip, Lp)>0 in Q- {t) 



ipt - F{D^ip, Dip, ip)>Om 0+{t). 



The above inequahties as well as (6.2 ) yield that tp is a, supersolution of (6.1 ) for < < < io 



if n is sufficiently large with respect to S. Since ip{x,0) = uo + 6 > uq, it follows from 
the comparison principle for ( |6.1[ ) that Un ^ <p for < t < to and for sufficiently large n. 
Therefore we have uj{x, 0) < uo{x) + S, and we conclude that uj{x, 0) < uq{x). 

To show that lo{x, 0) > uo{x) one can argue similarly as above to construct a subsolution 
of (6.1 1, with U in the proof of Lemma 4.1 and with the initial data uq — S. 



3. Now we construct n^ as follows. Note that, for any sequence rik -> oo, Proposition ro?2 
as well as the boundary data estimates given in steps 1.-2. yield that 

limsup* w„j. < u* . 

/c— ^oo 

For given Sk, choose Uk such that u„j_ satisfies 

Me,, - £k < Wnfc for < t < T, 



where u^^ is the minimal viscosity solutions of (1.1 1 with initial data Uq'': such nk exists 



due to Proposition 6.2 Now we conclude by applying Theorem 



4.3 



to u^ 



a 



6.4 Corollary. (For nonlinear F) Under the setting of Theorem 5.8. {u,i} converges to u 
in the sense of Corollary \6.3\ 



6.5 Corollary. Suppose uq £ V and F is as given in Coroll ary |5. 7| Then Un with the 
initial data Uq converges to the unique viscosity solution u of (1.1) with initial data uq in 
the sense of Corollary \6.. 



Appendix 

A Elliptic and parabolic theory 

In this section we collect a few classical results from the theory of fully nonlinear elliptic 
and parabolic operators used throughout this article. 



34 



A.l Proposition (Elliptic comparison). Let 17 C M" be a bounded domain and let F satisfy 
(1.6) and (1.7). Suppose that u € USC{n) is a viscosity subsolution and v S LS'C(f2) is a 
viscosity supersolution of the elliptic problem —F{D^u, Du, u) = in fl, and u < v on dil. 
Then u < v in Q. 

Proof. This is a consequence of ABP estimate applied to u — v, see (Aj Proposition 2.17]. D 

A. 2 Proposition (Parabolic comparison). Let S C M" xM. be a bounded parabolic domain 



and let F satisfy (1.6) and (1.7). Suppose that u e USC{Tj) is a viscosity subsolution and 
V G LSC(T,) is a viscosity supersolution of the parabolic problem Ut — F{D^u, Du, u) = on 
Y, such that u < v on dpY,. Then u <v inYi. 

Proof. See for instance |CIL Theorem 8.2] or |W, Corollary 3.15]. 



n 



A. 3 Proposition (Elliptic Hopf's lemma). Let u and v be respectively a subsolution and a 
supersolution of ~F{D^u, Du, u) — 0, u < v on a domain fl and u ^ v. If u ~ v at xq C^ dfl 
and Q has a interior ball touching xq, then 



y . f uixo) - u{xo - hC) . ,v{xo) -v(xQ-h^) 
lim mt ; > hm mi ; . 

/!,\,0+ h h\Q+ h 



(A.l) 



for any ^, ^- v > 0, where v is the unit outer normal to the interior ball at xq. provided that 
the right-hand side is finite. 



Proof. 1. For uj = u — v, we first show that 



liminf ^("°^-"J^"'^^^ >0. 
h\0+ h 



(A.2) 



Indeed, (1.6) yields 



-M- (M + N)-Si\p + q\- So \z + w\< F{-N, -q, -w) - F{M,p, z). 

Since —v is a subsolution of —G{D^v, Dv, v) — 0, where 

G{N,q,w)^~F{^N,-q,-w), 

we can use [AJ Proposition 2.14] to conclude that w = u + (— u) is a subsolution of 

-M^iD^u) - 6i \Dlj\ - So \uj\ = 0. 

Thus Hopf's lemma (AJ Proposition 2.15] applied to uj at xo (recall that a;(a::o) = 0) yields 
(A3. 



2. Using the property of liminf, we rewrite 

^^.^^ u{xo) ^ ujxQ - hi) ^ ^.^ .^^ v{xo) ~ vjxo - hi) 

h\0+ h h\Q+ h 

+ liminf "("°)'f"°-^^\ 

h\0+ h 



which together with (A.2) implies (A.l) 



a 
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A. 4 Proposition (Parabolic Harnack's inequality). Suppose that F satisfies the structural 
condition (1.6), is proper (1.7), and i^(0,0,0) — 0. Let ti, t2 satisfy < ti < t2. There 
exists a constant c, depending only on X, A, Si, Sq, ti,t2 and n, such that 

sup u < c inf u 

for any e G (0, 1) and any nonnegative continuous solution u of the parabolic problem on 
Qe = B,x (0,e2t2]. 

Proof. The main goal is to show that c does not depend on e. Let e G (0, 1) and let u > 
be a continuous solution of the parabolic problem on Q^ . Let us define the rescaled function 
u'^{x,t) = u{ex,e'^t) defined on Qi. We will show that we can apply Harnack's inequality 
to u^. 

Let iy9 be a smooth function and suppose that u"^ — ip has a strict maximum zero at 



{xo,to) G Qi. Then u — (p^'^ has a strict maximum at {exo,e'^to), and by (1.6) and (1.7) at 
(2^0,^0), 

£-2 {ipt - M+{D\) - Si \Dip\) < e-^ipt ~ e-^M+{D\) - e'^Si \D^\ 



Vt 



^ l/£ 



\l' -M^(p\^l')-Si 



Dtp 



lie 



F{D^ip^'',Dp^/',p^/') <0. 



Therefore u G S{Si, 0, 0) as defined in W 



Similarly, if w — 1^9 has a strict minimum zero at (xq, to)j then 
£-2 (<^^ „ M-{D^^) + Si \D^\ + <5o 1^1) 

> e-^ipt - e-^M-{D^p) + e^Mi \Dp\ + <5o \p\ 

i.e u G S{Si,So,0). 

Therefore we can apply the Harnack's inequality IWI Theorem 4.18] to w^, since its proof 



relies on the weak Harnack's inequality W Corollary 4.14] and the local maximum estimate 
W Theorem 4.16], and the function u'^ satisfies the hypotheses of both. D 



B Construction of barriers 



In this section wc construct necessary radially symmetric barriers for (1.1) used throughout 



the paper, first for fully nonlinear operators and then for the divergence-type operators. 



B.l Fully nonlinear operator 

Given a radially symmetric function p of the form 

ip{x) =■0(^1), 

its Hessian can be expressed as 

1 



^\ \ \x\ / a; 



36 



where {x (8) x)ij = XiXj. 
Note that 



{D^ip{x)){y) = j^tl;'{\x\)y for all y,y ±x 

Therefore the eigenvalues of D'^if{x) are 

ei-...-e„_i-^^'(|x|), e„ = V"(|a;|), (B.l) 

\x\ 

and the eigenvector for e„ is #t. 

B.l Example. Indeed, we can express the laplacian of a radially symmetric function as 

B.2 Example. If 7/;(p) is increasing and concave then the Pucci extremal operators can be expressed 



\x\ 
M + {D\) = Kin - l)^4l^ + V(N|). 



Of course, if V'(p) is decreasing and convex, we simply swap M^ and M . 

Now we state two propositions that guarantee the existence of barriers used in the proof 
of the comparison theorem in section [3] 



B.3 Proposition (Radial solutions). Let F in (|1.1[) satisfy (1.6 1 with given constants X, 



A, Si, do, and ( |1.8[ ). There is a constant p^ = pdX, A,Si,n) such that for any po G (Oj/oji 
d>0>6, d + 6>0, and uj > 0, there exists e > and a radially symmetric function 
(p{x,t), decreasing in \x\, and the following holds: 



(i) flip (resp. —fJ.(p) is a classical subsolution (resp. supersolution) of (1.1) for any /i > 
on 

K = {x : po- e <\x\ < po+e} X {-e, e), 

(ii) \Dip^\ = a and \Dip~\ = —b at \x\ = pa, t = 0, 
{Hi) the zero set moves with velocity Co, i.e. 

{x : (p{x, t) = 0} = {x : \x\ = po + tot} for t G (—£,£). 

Proof. To construct an increasing subsolution of the parabolic problem, we consider the 
radially symmetric positive function 

^{x, t)=ct + ^(p) -.^ct + a {p-'' - po^) + /3 (p2 _ pI) 
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for some a, /3, c, 7 positive, a and /3 are chosen so that ip is decreasing in p. A straightforward 
calculation yields 



^" = £,7(7 + l)p-T-2 + 2/3 > 0, 



1/,' = p (-a7p-'^-2 + 2;3) < 



Then using (B.l I and ( 1.5 1, we have at p = po f^^^d t = 0, 

ifit - F{D^ip,Dif,ip) <c-M-{D^ip) + S, \Dip\ i+So \^\ - F(0,0,0)) 
= c - a(A(7 + 1) - (n - 1)A - (5ipo)7P J~' 

-2/3(A + (n-l)A-<5ipo) 
= : c — ari — /3r2. 

The term r2 is positive for any po that satisfy 



(B.2) 



A + (n-l)A 
< Po < ^ =: Pc- 



(B.3) 



Furthermore, for po in this range, one can choose 7 large enough depending only on A, A, Si 
and n so that ti is also positive. If c < j3t2 as well, we observe that due to the continuity, ip 
with such parameters is a strict subsolution of the parabolic problem (and also the elliptic 
problem) in a neighborhood of {(a;,0) : \x\ = po}. 



Given a < 0, a; > and po satisfying (B.3 1, we can choose parameters a 
a way that \Dip\ 



and c in such 
Indeed, since ip is 



|a|. Dip ■ V = a and V^ = 10 ai t = Q and |a;| — po 
smooth, decreasing in p and increasing in i, we can express the normal velocity Vv on the 
zero level set {|a:| = po} at t = as 



K 



\D^\ 



2/3po - a7Po ^ 



The conditions V^ = ijJ and l-Diys] = |a| yield c = oj |a|. Then we find (3 large enough so that 
c < /3r2 and a < 2/3po and finally we solve for a > from \a\ = \Dip\ = 2/3po — a7Po ''^~ ■ 
Therefore (B.2) is strictly negative on {(a;,0) : |a;| = po}, and continuity yields that <p is 
in fact a strict subsolution of the parabolic problem on a neighborhood of this set. In 
fact, due to smoothness of ■0, it is possible to replace the term ct by the unique smooth 
increasing function of f, T{t)^ with t'(0) ~ c, which guarantees that {x : ip{x,t) = 0} = 
{x : \x\ ^ p + tut} for i in a neighborhood of i = 0, while at the same time (f is still a strict 
supersolution. Finally, since the right-hand side in (B.2 1 applied to fnp is 1-homogeneous in 



p > 0, pip is also a subsolution on the same neighborhood. 

The above construction also provides a subsolution of the elliptic problem for any c e M 
and /3 > 0. Moreover, function ip = —ip is a decreasing supersolution of the parabolic 
problem or a supersolution of the elliptic problem since A4^{M) = — Al+(— Af). D 



The free boundary of a solution of ( 1.1 1 can propagate arbitrary fast in some situations. 



In particular, if a supersolution is positive on a boundery of an open set G, then it will 



immediately become positive in G, as observed in the proof of Lemma 3.12 We show this 
using a barrier constructed in the following lemma in the form of the fundamental solution 
for the heat equation in one dimension.. 
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B.4 Lemma. Suppose that F satisfies ( |1.6[ ) with constants \,A,5i,So. Let G be a bounded 
open set and let c,S>Obe given positive constants. Then there exists a classical strict 
subsolution Lp of (1.1) on E, S := G x (0,(5], such that ip < c on I], (p{-,0) < on G and 
ip{-:S) > Q onG. 

Proof. We can assume that G 7^ 0. Let us define d := diam G. By translating G if necessary, 
we can also assume that G <Z {x : d < xi < 2d]. We shall consider (j) of the form 



{x,t) = il^{xi,t) -e, 



V'(a;i,t):=t-i/2 



exp 



Akt 



for suitable constants k,e > chosen below. 

First note that "(pi'it) is convex for x\ > 2kt. Then (1.6) and (1.5) yield 



0t - F{D^4>, D(f>, (j))<(j)t- M-iD^^) + 61 |D0| + So \^\ 



' {ipt - Xil^xixi + Si l^pxi \+So\ip- e\) 



< a 



{x\ - 2kt){k - \) Si\xi\ 



AkH"^ 



2kt 



So 



V'. 



By choosing fc, < fc <C min(^, A), we can guarantee that the bracketed quantity is negative 
for all te (0,2<5), xi G [d,2d]. 

Observe that ip{-,t) is decreasing when xi > 0. Moreover, for the fixed fc, there are 
77 £ (0, S) small and £ > such that 



-1/2 
7] ' exp 



-72 1 



Akrj 



<£<((5 + ry)-i/2 



exp 



k(S + ri) 



since the left-hand side goes to and the right-hand side is bounded from below as 77 — )■ 
0. With this choice of k,rj and e we see that 4>{-,rf) < in G C {x : d < Xi < 2d} and 
(t>{-,S + 'n) > on G. 

Now it is easy to verify that the function 



a 



ipix,t) = a0(x,t + 77)+ - -0(x,t + ?])_, 

where a > is sufficiently small so that ip < c on S, has the properties asserted in the 
statement of the lemma. D 



B.2 Divergence- form operator 

Next we consider the operator 



5(u)i-V-(*(6(w))Vu) =0 



(B.4) 



where b and ^ satisfy the assumptions from the introduction. Since this operator is not 
positively homogeneous of degree one in u as the Pucci operators, it is necessary to construct 



barriers analogous to Proposition B.3 in two steps. In particular, the barrier of Proposi- 
tion |B.5| is used in the proof of Lemma |3.21[ We do not present a modified version of 
Lemma B.4 since the barrier for problem (B.4) can be constructed in a similar way. 
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B.5 Proposition. For given lo > 0, po > and M > there exists rjQ > such that for all 



rj € (0,770) there is a strict classical supersolution ip of the divergence form problem (B.4) 



E := {(x, t) : po+oot < \x\ < po + ojt + rj, t e (— t, 00)} 

with T = ^ if LO > Q and t — 00 if uj — 0. Furthermore Lp{x,t) = for \x\ = po + tot and 
ip{x, t) > 2M for \x\ = po+UJt + rj. 



Proof. Let us set 



ujb'(s) 2(n-l) 

ki := max ^ ,, , ,\ + ^ > 0, 

se[o,3M] *(6(s)) Po 



kn ;= max 



W{b{s))\b'{s) 



se[o,3M\ ^{b{s)) 



(B.5) 



<2M. 



(B.6) 



Set rjo := k^ . Now we can fix any 77 e (0, 770), and for such 77 choose k > k2 large enough 
so that 

k-k2 1 1,-1, k-k2 

-fcr--r>^ and k log-r— 

Finally, we choose a > 1 large enough such that 

2Af < k'^ log(a/c77 + 1) < 3Af. 



Note that this is always possible since \jr^ — -^ ^^ Vo > V ^^^'^ ^ ^ log ^'^-^ 
and fc77 + 1 < ^^. 



fci 



We will show that the function ip of the form 

(p{x,t) =ij{\x\-ujt- Po), ip{s) 



log(afcs + 1) 
k ' 



as fe — )■ c», 



(B.7) 



is the sought supersolution of (B.4) 



Indeed, recall that |a;| > po/2 and ip{x, t) G [0, 3Af ) in S. With the help of fci, fc2 defined 
in (B.5), we can estimate 



*(%)) 



(%)t-V-(*(%))V^)) 



ujb'{^) ,, *'(5(^))6'(^) 



*(%)) 



V- 



*(K¥')) 



v/^- 



71- 1 



^' - iP" 



> -V''(fcl + k2^p') - V". 



A straightforward differentiation of tp defined in (B.7) then yields 

-i;\s){k, + k2tP\s))-^"{s) 



ki +k2 



aks + 1 \ aks + 1 

a f a{k — k2) 



a^k 



{aks + 1) 



aks + 1 V ciks + 1 



fci > 0, 



for all s = \x\ — tot — Po e (0,77) due to the choice of k in (B.6) and a > I 



n 
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The barrier constructed in Proposition B.5 is suitable for the construction of a barrier 



for the Richards problem similar to Proposition B.3 



B.6 Proposition. For any po > 0, a > > b, a + b > 0, and a) > there exists e > and 
a radially symmetric function ip{x,t), decreasing in \x\, and the following holds: 



{i) (p is a strict classical subsolution of (1.1) in the divergence form (B.4) on 

K ^{x : po~e <\x\ < po+e} X (-e, e), 

(ii) \Dip\ ~ a and \Dip\ = —b on \x\ = pa + ujt, 
{Hi) the zero set of tp moves with velocity to, i.e. 

{x : ip{x, t) —Q} = {x : \x\ = Po + <^t}- 
Proof. We will use the function constructed in the proof of Proposition |B3] Let w = a) and 



M — 1. Choose any k > ^2, where fc2 was defined in (B.5), and set a = a. Let -0(5) be as 
defined in (B.7) and for some ^ > and t g {—S,S) let us set 



(p{x,t) 



(2-n)(po+ioty- 



for s e (-(5,0], 
s^ for s > 0, 



where s = \x\ — po — ut. Then due to the computation in that proof, we can find (5 > such 
that 

• '0('S) is well-defined and — V'('S) < 1 = ilf for s E {—S, 0), 



(p is a strict classical subsolution of (B.4) in 



{{x, t) : \x\ -po-ute i-5, 0), t e (-<5, <5)}, 

• [IJi^^l = a, [IJi^^l = —b on |a;| = po + tot. 
We finish the proof by choosing e G (0, 5) small enough so that 

K C {(x, t) : |a:| - p^ - ijjt e (-(5, oo), t e (-(5, 5)}. 
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